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It is shown that the norm corrections, introduced to avoid the violation of the constraints 
on the depletion of the hole states in the standard perturbative 2p2h approach, leads 
in nuclear matter to a dependence of the momentum distribution with the total nucleon 
number. This unphysical behavior, which in turn makes the depletion to be non-extensive, 
arises from contributions of disconnected diagrams contained in the norm. It is found 
that the extensivity is again recovered when the ApAh excitations in the ground state are 
included, and a reasonable value for the total number of nucleons promoted above the 
Fermi level is obtained. 
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I INTRODUCTION 



A great amount of theoretical effort was devoted to study the influence of ground-state 
correlations on the nuclear physical observables. In most of them the nuclear response 
function is evaluated after introducing two particle - two hole {2p2h) admixtures within 
first-order perturbation into the OpOh ground state wave function |T|, |, ^ ^ Q . It was 
found that 2p2h ground-state contribution modify significantly the mean field results for 
the strength function. Recently, however, the use of this perturbation procedure has been 
objected by Van Neck et al, [0 because it largely overestimates the effect of ground- 
state correlation, when the norm corrections are neglected. This happens because the 
perturbation adds a very large number of relatively small excited 2p2h components, with 
relative weights in the perturbed wavefunction typically large enough numerically as to 
strongly enhance the norm corrections. Specifically, it is reminded in Ref. that a 
shell-model approach for occupation numbers in nuclei imposes constrains in the sense 
that the number of particles lifted out the Fermi sea has to be < 2. This constrain is not 
respected in the above mentioned works as put in evidence by several calculations of the 
occupation numbers within the same perturbation approach |]10|. 

In a recent study of the Gamow- Teller strength the ground-state was worked out 
within a finite nuclei formalism up to the first order, both in the Rayleigh-Schrodinger 
and Brioullin-Wigner schemes, and including the effects of the norm as was suggested 
in Ref. 0. The result was the strong reduction of 2p2h contribution to the response 
function. In the present work we intend to implement better controlled approximations 
for the correlation structure of the ground-state for the case of nuclear matter, which may 
have important consequences in the evaluation of the response function in the quasielastic 
region. As a first step in this direction we analyze here the momentum distribution in the 
ground state. The reason for that is that, as shown within a simple model by Takayanagi 
0, the longitudinal response in the quasi free region is directly related to the momentum 
distribution of nucleons in nuclear matter. Our initial intention was to include the norm 
effects in the some way as we have previously done for the finite nuclei ^j. But, we 
soon discovered that this procedure makes the depletion of the nuclear core not to be 
anymore an extensive quantity because of the momentum distribution dependence on 
the total nucleon number A. Clearly, such a behavior is unphysical and has its origin 
in the contributions arising from the infinite series of disconnected Goldstone diagrams 
contained in the norm, as can be put in evidence by performing a perturbation expansion 
of the later. In order to exactly cancel these disconnected graphs one is forced to enlarge 
the configuration space by including npnh amplitudes with n > 2. This is done in the 
present work for the 4p4h components in the ground state wave function. To go beyond 
this configuration space is not only a cumbersome task but hard to justify also. 
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II FORMALISM 



The occupation number for the single-particle state k in the ground state |0) is defined 

as 

n{K) = {0\J(K)a(K)\0), (2.1) 

and the total number of particles promoted above the Fermi level, A^>, and those remaining 
below the Fermi level, A^<, are 

N< ^ n{K) ^A-N>. (2.2) 

Within the Hartree-Fock (HF) approximation, |0) = \HF) = \OpOh) one obtains the well 
known step function for n(/t), i.e., n^n) = 6{tp — e^), where and are, respectively, 
the single particle and Fermi energies. This is the zero order approximation for the 
occupation numbers distribution. When npnh correlations are added to the ground state 
wave function, states above the Fermi level are populated with the corresponding depletion 
of the nuclear core, and the occupation number takes the form 

n{K) = e{e^ - e J + 5n{K) , (2.3) 

being 

7V>= ^ 5n{K). (2.4) 



A 2p2h correlations 

The ground state wave function, with the first order 2p2h perturbations included to 
the HF state, reads 



p's,h's 



(2.5) 



where 



(j)iP2hMV\HF) 



'PiP2hih2 



(2.6) 



2 



are the 2p2h amplitudes, and 



Ar=^_=, (2.7) 

Y J- + J^2p2h 



is the overall normalization factor, with 

J^2p2h = - \(^PlP2hih2f ■ (2-8) 

p's,h's 

\p1P2h1h2) and Ep-^p^hih2 stand, respectively, for the 2p2h configurations and single-particle 
energies, and V is the residual interaction. 

The usual second-order approximation for n{K,) is obtained by retaining only the first 
term in the expansion 

^f' = l-^f2p2h + ^^2p2h-■■■, (2.9) 

and for 6n{K,) one gets 



^ p's,h's \_^PiP2hih2 



V 2 

^PiP2hih2 



(2.10) 



where ^ip2/ii/i2 stand for the antisymmetrized matrix elements of V . The contribution of 
5n(K)*^^) is schematically shown in Fig. |I|, and 

= 2U2p2h. (2.11) 

This result is only valid when M2p2h 1, condition which is not fulfilled in most of the 
numerical calculations performed so far. This was precisely the reason why Van Neck et 
ai, [0 have proposed to use the second order normalized approximation (with A/" given 
by (p.7|)). This leads to the variation in the occupation numbers 



6n^^^\K) = Af^5n^^\K) (2.12) 
and the corresponding number of promoted particles is 

iVr^ = ^f^Ni'^ = 2 -^'f; < 2, (2.13) 

1 + M2p2h 

which is a satisfactory bound for A^>. Yet, StS'^^^k) contains an infinite series of 
disconnected diagrams, illustrated in Fig. |^. We show below how they are removed 
by enlarging the configuration space. 
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B 4]?4/i correlations 

The ground-state wavefunction with the 2p2h and 4p4/i admixtures included is written 



as 



|0) = M 



\HF) + 



1 



('21)2 XI Cp^P2hih2\PlP2hlh2) 
^ ' p's,h's 



+ 



\ p's,h's 



Cpip2P3Pihih2hzhi \PlP2P3P4hlh2h3h4) 



where 



^PlP2P3P4hlh2h3h4 



and 



{HF I V\piP2hi h2 ) {piP2hi h2 \ V\piP2PzPihi ^2 ^3 ^4) 



1 + ^2p2h + J^Apih 



with 



4p4/i 



E 



-PlP2P3P4'l.l/i2fc3'l4 I ) 



('41)2 ^ 
V '/ p's,h's 

which is the norm factor coming from the ApAh contributions. 
The occupation numbers in the ground state (|2.14|) read 



6n{K) = 



where 



(A\\2 E ("^K.Pi 
V ■/ p's,h's 



p-iP2PzP4hih2hzhi I 



Now, making an expansion up to the fourth order 

5n{K) = [1-N2p2h] + 5n^''\K) 

and after working out the ApAh amphtudes we get 



(2.14) 



(2.15) 



:2.16) 



(2.17) 



(2.18) 



:2.19) 



(2.20) 



(2.21) 
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or 



(2.22) 



with 



X 



p's,h's ^PUnhihi^paPAh-ihi 

(^^K,pi ~l~ ^K,P2 ^K,hl ^K,h2 



^PlP2hlh3^pZP4,h 2hA ^piP3hlh2^P2P4h3h. 



E. 



piP2hlh3Ep^p^h^fl4 



^PlPshlh2Ep2P4h3h4 _ 



I ■ PlPahihs ^P2P4h2h4 1 ^pp2h'ih4^p3P4hih2 



E. 



PlPihihs Ep2P4h2h4 



2 -En 



pP2h3h4^P3PAhih2 



(2.23) 



being the contribution of the connected diagrams illustrated in Fig. ^ The diagrams (a), 
(b) and (c) describe the different contributions that can arise from the first term, while 
the graph (d) represents the second term. For the depletion number we now get 



2M 



2p2h 



ipih 



2M 



2p2h 



[2.2A) 



where the last term comes from 5rS'^^\K). Note that the result (|2.24]) corresponds to the 
expansion up to fourth order of the quantity 



'2M2p2h + 



ipih 



1 + J^2p2h + AT, 



< 4 



' ipAh 



(2.25) 



which is the total depletion number obtained from the exact occupation probabilities 
( p.l8| ), and same as (|2.13|) , is also properly bounded. 

Regarding the results of this subsection, it should be remarked that: 

(a) In the expansion (p.9|) are retained the first two terms. 

(b) The contributions to 6n{K,) coming from the disconnected graphs brought around by 
J^2p2h are not present in final result because they cancel out with topological similar 
contributions arising in the wave function (|2.14|) , i.e., the first term of Sti^^^k) in 
Eq.(^). 

(c) The norm term M^pih does not contribute in a calculation up to the fourth order 
and has been ignored in the expansion of N"^. 

(d) We do not include the coupling among the 2p2h states ((2p2/i|\/|2p2/;,') = 0), because 
it would only lead to a redistribution of the 2p2h occupation probabilities, which 
does not have any effect on 5rS^\k). 



5 



C Evaluation of the occupation numbers distribution in 

nuclear matter 

In nuclear matter the HF ground state is approximated by the Fermi gas, the single 
particle quantum numbers k are the momentum k, the spin projection and the isospin 
projection mt of the particle, and n(/t) turns into the momentum distribution 

n{k) = e{l-k) + 5n{k), (2.26) 

where k is measured in units of the wave number kp that defines the Fermi surface. The 
corresponding depletion number now reads 

dkk^e{k-l)Sn{k). (2.27) 
The residual interaction is generically expressed in the form p 

nq) = E^'(?)oi(q)-02(-q) (2.28) 

where V^{q) are the interaction strengths and I = T, S, J encompasses the isospin, spin 
and total angular momentum quantum numbers of the operators 0^(q), defined as 

0°°°(q) = 1; 0°i°(q)=^(q-^); 0''\q) = x a) , 

Oio°(q) = r; 0"°(q) = ^(q ■ cr)r 0"^(q) = (q x cr)T. (2.29) 

The exchange contributions to 6n{k) will be dropped out since, as pointed out Van Order 
and Donnelly ITlI], they are small in comparison with the direct ones. In this way we get 



6n^^\k) = dq[9{k- 1)9(1 -\k + q\)J^(k ■q + q^,q) 

+ 6(1 - k) ^(|k - q| - l)^(k ■ q, q)] Y.{v\q)f{2T + 1)2^, (2.30) 

for the second order correction to n{k), 

M2p2h = —] dMq9{k - l)e{l - |k + q|)^(k ■ q + q) Y,[v\q)]\2T + 1)2-', (2.31) 

for second order correction to the norm, and 

5n^^^\k) = -^Jd\JdpJdqe{l-\l + q\)e{t-l) 

X [e{k-i)g{k,\,p,q) + e{i~k)g{k-q,i,p,q)], (2.32) 
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for the fourth order correction to n{k). In the above equations 



2P 



(2.33) 



9{x-l)9{l 



q|) 



+ 



2tt 

q 

In 



(a + q ■ x)^ 
-1 + ln2 + ln 



(ix 



q — a + q 



q' 



2a 



a 

- — In 



q — a + q 



q — a + q 



q^ — a — q 



aq— a + q^2 
q^ q^ — a — q q 



q — a 
%-2)|, 



d{2 - q) 



[2.M) 



and 



+ 



g{\i, 1, p, q) = [^(1 - |i + p|)^(i ■ p + /2, /) + e{i - |k + p|)^(k ■ p + /) 
^(|i + q + p|-im-i-p-q-p,0 + ^(|k + q + p|-im-k-p-q-p,/)] 

^^(l-|k + q|) 



(k ■ q + g2 + 1 . q) 



5:b^(g)]2(2T + l)2^5:K(/)]2(2r + 1)2^'. 



2.35) 



In Eq. (|2.32|) we do not include the contributions arising from the second term in Eq. 
( p.23| ), because, as discussed later on, they are relatively small. 

From Eqs. (|2.11| ), ( p.l3|) and (|2.31|) it can be easily seen that A^>^'' fulfills the exten- 



sively condition, but N^^\k) does not. This unphysical behavior of N^"' is due to the 
disconnected diagrams contained in dnP'^^k) and is circumvented by including the ApAh 
excitations. This lead us to the result in Eq. (|2.24|) , which again satisfies the requirement 
on extensivity. 



{2N) 



III NUMERICAL RESULTS AND 
CONCLUSIONS 

The momentum distribution n{k) has been evaluated for the normal nuclear matter den- 



sity, i.e., kp = 1.36 fm, and by parameterizing the residual interaction (|2.28|) as follows: 



V'''{q)=C^{q)f, V^'%q) = CMf\ V''\q) = V''\q) = C^{q)g 



V'''{q) = CM 9' 



q^ 

q^ + ml 



V''\q)=CMg', 



(3.1) 
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Here /, /', g and g' are the Landau-Migdal parameters, C.„{q) = f^/^lTl{q), = ^^/kp 
and 

/\2 _ ™2 

r.(^) = (3-2) 

is the ttNN form factor. 

We do not incorporate exphcitly the p meson since, as pointed out by Dickhoff |T^, 
it produces too much suppression of the tensor force. Yet, the contribution of heavier 
mesons is taken into account empirically by fixing the values of g' and (as done by 
Alberico et al, P]) to reproduce the effective tensor and central components of the force 



derived by Dickhoff from a G-matrix calculation (|T^, |T3]). This procedure yields g' = 0.5 



and Att = 800 MeV/{hckp). On the other hand, for /, /' and g we have found more 
appropriate to adopt the values obtained by Backman et. al. from the analysis of 



stability conditions and empirical values of physical observables, i.e., f = —0.12, /' = 0.39 
and g = 0. 

The following three different ways of adding the ground state correlations to the Fermi 
gas distribution will be numerically compared: 

1. The usual 2p2h second order approximation, i.e., Eq. ( 2.30| ) is employed in the 
evaluation of 6n{k). 

2. The above approach is modified by including the effect of the norm through Eqs. 
( pJ2|) , Q and 



3. Fourth order approximation within the 2p2h + ApAh subspace, in which 5n{k) is 
given by Eqs. and (^L^ ). 



The momentum space integrals have been evaluated for energies up to 300 MeV . 
This limit has been established by examining the behavior of N^/A (with Sn{k) given by 
( p.30| )) as function of the cutoff energy. The results are displayed in Fig. |[ 

The resulting momentum distributions n{k) are shown in Fig. |^. In the approximation 
2), where the results are A dependent, we have chosen, as an example, the ^°Ca nucleus. 
For this case the empirical occupation numbers are available from an optical potential 
analysis of experimental cross sections (see also ref. 0). The second order approach 
produces a pronounced high-momentum tail above the Fermi level with a corresponding 
strong reduction of the low-momentum part as compared with the Fermi gas. In the 
normalized case 2), the previous variation of momentum distribution 6n{k) is just rescaled 
by the factor A/"^ = 0.11 and this diminishes drastically the effect of the ground state 
correlations. (Note that for A = 90 N"^ = 0.05.) Within the approximation 3) the 
depletion of the momentum distribution is quite significant, but still relatively small in 
comparison with that obtained in the plain 2p2h approach. This effect can be interpreted 
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as the reduction of the 2p2h occupation probabihties by the couphng with ApAh states. 
It is worth noting that the effect of the last term in Eq. ( p.23| ) cannot be distinguished 
visually in Fig. ^ 

The calculated numbers of particles shifted above the Fermi level, for the implemented 
approximations, turn out to be A^> = 16.8, 1.8 and 1.9, respectively. The last result is, 
not only physically sound, but also consistent with the above mentioned study [1^, which 
yields A^>(^°Ca) = 2.7. 

We summarize our conclusions as follows: 

• The plain second-order approximation for the momentum distribution grossly over- 
estimates the effect of the 2p2h ground-state correlation. This is put in evidence by 
the very large result obtained for the number of particles shifted above the Fermi 
level, as already pointed out 

• The introduction of norm corrections, following the recipe given in Ref. 0, leads 
unavoidably to an A dependent momentum distribution, which is engendered by 
the disconnected diagrams present in the norm. As a consequence the number 
depletion becomes a non-extensive quantity for the nuclear matter, which is very 
strong drawback of the approximation. 

• We have found out that the just mention unphysical contributions of the discon- 
nected diagrams are cancelled by the inclusion of 4p4h correlations. Besides, the 
contribution of the ApAh connected diagrams strongly hinder the effect of 2p2h 
correlations, yielding a momentum distribution with properly bounded depletion 
number. 

• The present results seem to indicate that the 2p2h and ApAh correlations are the 
dominant degrees of freedom to be considered in the description of the nuclear 
ground state. In addition, we feel that the evaluation of the momentum distribution 
up to fourth order is a reasonably good approximation and that is not crucial the 
inclusion of higher order correlations beyond the ApAh ones. 



This work was supported in part by a Fundacion Antorchas grant. 
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Figure Captions 



Figure 1: Goldstonc diagrams corresponding to the second-order correction Sn^^'' (k). Each hne indicates 
schematically a particle or a hole state, the dots represent the residual interaction and the encircled dots 
correspond to the number operator. 

Figure 2: Goldstonc diagrams corresponding to the expansion of the norm in the normalized sec- 
ond-order 2p2h approximation 6n^'^^\K). 

Figure 3: Goldstonc diagrams corresponding to the different contributions in the fourth order 4p4/i 
correction 5n^^'-'\K). 

Figure 4: Plot of the second-order 2p2h approximation for Ny/A, as function of the cutoff energy Ec 
(measured in MeV) for 2p2h excitations. 

Figure 5: Momentum distribution n(«;) for the different approximations: Fermi gas (thin full lines), 

second-order approximation with 2p2h correlations {5n{K) = Su^'^^k)) (long-dashed lines), normalized 
second-order approximation for ^"Ca {Sn{K) = 5n'-^^^{K)) (short-dashed lines), and fourth-order approx- 
imation with 2p2h + 4p4/i correlations, {6n{K) = (5n^^'(K) -|- Sn''^'~^\K)) (thick full lines). 
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